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Abstract: The simplest unified extension of the Minimal Supersymmetric Standard 
Model with bilinear R-Parity violation provides a predictive scheme for neutrino masses 

and mixings which can account for the observed atmospheric and solar neutrino anomalies. ^ 

I Despite the smallness of neutrino masses R-parity violation is observable at present and [\3 

I future high-energy colliders, providing an unambiguous cross-check of the model. ^ — ^ 

o 

00 



1. Introduction 

The announcement of high statistics atmospheric neutrino data by the SuperKamiokande 
collaboration [1] has confirmed the deficit of muon neutrinos, especially at small zenith 
. angles, opening a new era in neutrino physics. Although in the past have been considered 

alternative solutions for the atmospheric neutrino anomaly [2, 3] it is now clear that the 
^ , simplest interpretation of the data is in terms of to Vr flavor oscillations with maximal 

^ I mixing. This excludes a large mixing among Ur and fg [1], in agreement also with the 

CHOOZ reactor data [4,5]. On the other hand, the persistent disagreement between solar 
neutrino data [6] and theoretical expectations [7] has been a long-standing problem in 
physics. Recent solar neutrino data from the Kamland collaboration [8] and the latest 
results from the SNO collaboration [9], clearly indicate that we have MSW conversions 
with a large mixing angle [10]. For the solar neutrino parameters we get for the best fit 
point [10]: 

tan^ 6lsoL = 0.43, Aml^^ = 6.9 x 10~^ eV^ , (1.1) 

confirming that the solar neutrino mixing angle is large, but significantly non-maximal. 
The 3(7 region for 9 is: 

0.30 < tan^ O^ol < 0.64, (1.2) 
while the 3a region for Am^Q^ range is given by, 

5.4 X 10"^ eV^ < Am^oL < 9.5 x 10"^ eV^ (1.3) 

* Speaker. 
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On the other hand, current atmospheric neutrino data require oscillations involving i/,- 
[1]. The most recent global analysis gives [10], 

sin^ 6'atm = 0.52 , Aml,^ = 2.6 x 10"^ eV^ (1.4) 

with the 3a ranges, 

0.31 < sin^ 6Iatm < 0.72 (1.5) 

1.4 X 10"^ eV^ < Am^^^j < 3.7 x 10~^ eV^ . (1.6) 

Many attempts have appeared in the literature to explain the data. Here we review 
recent results [11-16] obtained in a model [17] which is a simple extension of the Mini- 
mal Supersymmetric Standard Model (MSSM) with bilinear R-parity violation (BRpV). This 
model, despite being a minimal extension of the MSSM, can explain the solar and atmo- 
spheric neutrino data. Its most attractive feature is that it gives definite predictions for 
accelerator physics for the same range of parameters that explain the neutrino data. 



2. The Model 

Since BRpV has been discussed in the literature several times [11, 12, 17-20] we will repeat 
only the main features of the model here. We will follow the notation of [11, 12]. The 
simplest bilinear model (we call it the MSSM) is characterized by the superpotential 

W = Wmssm + Wi^^ (2.1) 

In this equation Wmssm is the ordinary superpotential of the MSSM, 



W = eab 



4QtU,Ht + h%Q\D^H2 + h%L\R,H2 - fiH^Ht (2.2 



where i,j = 1,2,3 are generation indices, a, 6 = 1,2 are SU{2) indices. We have three 
additional terms that break R-parity, 

W^^ = eabtiLIHl. (2.3) 

These bilinear terms, together with the corresponding terms in the soft supersymmetric 
(SUSY) breaking part of the Lagrangian, 

Vsoft = Kt/f " + eabB,e,LtHt (2.4) 

define the minimal model, which we will adopt throughout this paper. The appearance 
of the lepton number violating terms in Eq. (2.3) leads in general to non-zero vacuum 



expectation values for the scalar neutrinos (Pj), called Vi in the rest of this paper, in 
addition to the VEVs vu and vd the MSSM Higgs fields and H^. Together with the 
bilinear parameters the Vi induce mixing between various particles which in the MSSM are 
distinguished (only) by lepton number (or R-parity). Mixing between the neutrinos and 
the neutralinos of the MSSM generates a non-zero mass for one specific linear superposition 
of the three neutrino flavor states of the model at tree-level while 1-loop corrections provide 
mass for the remaining two neutrino states [11, 12]. 
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3. Neutrino Masses and Mixings 

3.1 Tree Level Neutral Fermion Mass Matrix 

In the basis = (— iA', —iX^, H^,H^, Ug, v^, Vr) the neutral fermions mass terms in the 
Lagrangian are given by 

'Cm = -^(V'°fMivV'° + /i-c. (3.1) 



where the neutralino/neutrino mass matrix is 

Mat : 



m 



(3.2) 



with 



Ml 






M2 







m 



ai 

0.2 

as 



(3.3) 



where = {—\g'vi, ^gvi,Q,ei). This neutralino/neutrino mass matrix is diagonalized by 



3.2 Approximate Diagonalization at Tree Level 



(3.4) 



If the parameters are small (we will show below that this is indeed the case), then we 
can block-diagonalize Mat approximately to the form diag(me//, A^^o) 



1 T + 

rueff = -m ■ M o m = 



4 det(A^ 0) 



AeAf, Al Af,Ar 
\ AeAr A^Ar A^ J 



with 



Aj = iivi + Vdei . 
The matrices A'' and V^, diagonalize ■M-)^o and me// 

N*M^oN^ = diag(m^o) ; yJme//K = diag(0, 0, m,,), 

where 



Tr{meff) = 



(3.5) 



(3.6) 



(3.7) 



(3.8) 



4 det(7W^o) 

So we get at tree level only one massive neutrino, the two other eigenstatcs remaining 
massless. The tree level value will give the atmospheric mass scale, while the other states 
will get mass at one loop level. Therefore the BRpV model produces a hierarchical mass 
spectrum. We will show below how the one loop masses are generated. 
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3.3 One Loop Neutrino Masses and Mixings 
3.3.1 Definition 

The Self-Energy for the neutrahno/neutrino is [12], 




[PL^f^ + Pr^^] - [PLllf- + PrU^] } 



Then the pole mass is, 
with 



AM, 



1 /'ttV/_2n , -nVr^2^\ 1 



A=0 



where 



- i (s^ + s^) ; = i (n^ + n«) 



an the parameter that appears in dimensional reduction is, 

2 



A 



A-d 



+ ln47r 



(3.9) 

(3.10) 

(3.11) 
(3.12) 

(3.13) 



3.3.2 Diagrams Contributing 

In a generic way the diagrams contributing are given in Fig. 0, where all the particle in the 
model circulate in the loops. These diagrams can be calculated in a straightforward way. 
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Figure 1: Diagrams contributing at one loop to the neutrino mass matrix 

For instance the W diagram in the ^ = 1 gauge gives 
1 ^ 



CD 

O 
O 
00 



fc=i 



2 2 2 \ 



(3.14) 



fc=i 



where Bq and i?i are the Passarino-Veltman functions [21], and O™^, O'^^^ are coupling 
matrices. Explicit expressions can be found in [12]. 

3.3.3 Gauge Invar iance 

When calculating the self-energies the question of gauge invariance arises. We have per- 
formed a careful calculation in an arbitrary gauge and showed [12] that the result was 
independent of the gauge parameter ^. 



- 4 - 



International Workshop on Astroparticle and High Energy Physics 



Jorge C. Romao 



4. Results for the Solar and Atmospheric Neutrinos 



4.1 The masses 



> 

(D 



The BRpV model produces a hi- 
erarchical mass spectrum for almost 
all choices of parameters. The largest 
mass can be estimated by the tree 
level value using Eq. (|3.8| ). Correct 
Am^TM can be easily obtained by an 
appropriate choice of |A|. The mass 
scale for the solar neutrinos is gen- 
erated at 1-loop level and therefore 
depends in a complicated way in the 
model parameters. We will see be- 
low how to get an approximate for- 
mula for the solar mass valid for most cases of interest. Here we just present in Fig. |2|, 
for illustration purposes, the plot of the three eigenstates as a function of the parameter 
e^/|A|, for a particular values of the SUSY parameters, mo = /U = 500 GeV, tan/3 = 5, 
B = —A = ruQ. For the R-parity parameters we took |A| = 0.16 GeV, 10 * Ag = = Al- 
and ei = £2 = €3. 




Figure 2: neutrino masses 
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4.2 The mixings 

Now we turn to the discussion of the mixing angles. As can be seen from Fig. |2[ if 
e^/|A| <C 100, then the l-loop corrections are not larger than the tree level results. In this 
case the flavor composition of the 3rd mass eigenstate is approximately given by 



Ua3 ~ A„/|A| 



(4.1) 



As the atmospheric and reactor neutrino data tell us that — > Vr oscillations are preferred 
over f/t — > fe, we conclude that 



Ae <C A„ ~ A^ 



(4.2) 



are required for BRpV to fit the data. This is sown in Fig. ^ a). We cannot set all the Aj 
equal, because in this case J7e3 would be too large contradicting the CHOOZ result as shown 
in Fig. ^ b). We have then two scenarios. In the first one, that we call the mSUGRA case, 
we have universal boundary conditions of the soft SUSY breaking terms. In this case we 
can show [11,12] that 

Cr Ar 



A„ 



(4.3) 



Then from Fig. |3|b) and the CHOOZ constraint on C/^g, we obtain that both ratios in Eq. (^ 
have to be small. Then from Fig. ^ we conclude that the only possibility is the small angle 
mixing solution for the solar neutrino problem. In the second scenario, which we call the 
MSSM case, we consider non-universal boundary conditions of the soft SUSY breaking terms. 
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Figure 3: a) Atmospheric angle as a function of |A^|/ y^A^ + A^. b) C/^g as a function of 
|A,|/^A2+A2. ^ 

f , > 

We have shown that even a very smah deviation from universahty (less then 1%) of the 
soft parameters at the GUT scale relaxes this constraint. In this case 

A, 



2 



(4.4) 



Then we can have at the same time small U^^ determined by Ag/A^ as in Fig. ^ b) and 
large tan^(0soL) determined by eg/e^ as in Fig. ^ b). After the Kamland and SNO salt 
results, this is the only scenario consistent with the data. 
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Figure 4: Solar angle as function of: a) |Ae|/yA^ + A^ ; b) ee/c/i- 
5. Approximate formulas for the solar mass and mixing 

In all the previous analysis we used a numerical program to evaluate the one-loop masses 
and mixings. It is however desirable to have analytical approximate results that can give 
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us quickly the most important contributions. We have identified that these are the bottom- 
sbottom loops and the charged scalar-charged fermion loops. Then, by expanding in powers 
of the small R-Parity breaking parameters ej, we get approximate formulas for the solar 
mass scale and mixing angle as we will explain below. 

5.1 Bottom-sbottom loops 

The contribution from the bottom-sbottom loop can be expressed as [14], 



AM,-. 



NcTHb, 



eiej 



(5.1) 



where are the e, in the basis where the tree level neutrino mass matrix is diagonal, 



e,; 



the a, are functions of the SUSY parameters, and 

^bI^''' = Bo{0, ml ml) - Bo{0, ml ml) . 



(5.2) 



(5.3) 



The different contributions can be understood as coming from different types of insertions 
as shown in Fig. |^. In this figure open circles correspond to small R-parity violating 
projections, full circles to R-parity conserving projections, and open circles with a cross 
inside to mass insertions which flip chirality. With this understanding one can make a one 
to one correspondence between Eq. (|5.lD and Fig. |5|. 





a3|A|(5j3 + 6ej b az\A.\5iz + hit 



Figure 5: Bottom-sbottom loop and different types of insertions. 



5.2 Charged Scalar-Charged Fermion Loops ^ 

Here the situation is much more complicated. As the charged scalars are now combination^^ 
of the six sleptons and of the two MSSM charged Higgs bosons we have many more particles 
that can be exchanged. For instance, the contribution from the staus is given in Fig. ^. We 
see that, even for the staus, we have more diagrams than for the case of the bottom-sbottom 
loop. This is due to the fact that the breaking of R-parity is in the leptonic sector of the 
theory. So we can have R-parity violating insertions both on the external neutrino legs 
as well as in the particles circulating in the loop. A complete description of all the terms 
contributing to this loop can be found in Ref. [14], where one can also find the analytical 
expressions. 
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Figure 6: Charged scalar loops: the stau contribution 



5.3 Analytical vs Numerical results 



O 



We now compare the approximate analytical formulas with the full numerical calculation. 
This is shown in Fig. ^. On the left panel we consider a data sample where the neutralino 
is the Lightest Supersymmetric Particle (LSP). Shown are also the bands for LOW and ^ 
LMA solutions of the solar neutrino problem, as well as the much narrower band obtained > 
after the recent Kamland [8] and SNO salt [9] results. On the right panel we have the same > 
situation for a data sample where the scalar tau is the LSP. We see that the agreement ) 
is specially good (better than 10%), just in the region which is compatible with the most 
recent data. 
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Figure 7: Comparison between the analytical and numerical results for the solar mass. In the left 
panel the LSP=x*' while in the right panel the LSP=f. The red (dark) band corresponds to the 
latest neutrino data [10]. 
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5.4 Simplified approximation formulas 

In the basis where the tree-level neutrino mass matrix is diagonal the mass matrix at 
one-loop level can be written as 



'^cieiei cieie2 cieies ^ 

Ci?2?l Ci?2?2 Ci?2?3 (^-^^ 

\c1e3e1 Ci?3?2 Co|Ap + Ci?3?3 J 



where 



= 4det(>l,o) ^'-'^ 
ci = sin(2^^) mb AbI^'^ (5.6) 

The dots in Eq. ( |5.4[) correspond to the terms that are not proportional to the e, x ej struc- 
ture, as can be seen from Eq. Assuming that the bottom-sbottom loop dominates, 
the €i X ij structure is dominant, and the matrix can be diagonalized approximately under 
the condition 

X ^ ^ « 1 (5.7) 
co|A|2 

Then we get 

m,, ^ sm{2e~,)m,ABl^^' ^^^^ (5.8) 

and 

tan2esoL = 4 (5-9) 



^2 

The results for the masses are presented in Fig. ^. As in Fig. ^, on the left panel 
we have the data set where the neutralino is the LSP while on the right panel the LSP is 
the stau. We see that the agreement is fairly good, particularly in the region allowed by 
the present data. For the solar angle the results are shown in Fig. |9|. We see that the 
agreement is not as good as for the masses, even if we restrict ourselves to the present 
allowed values of tan^ 6'soL 1 shown by the red (dark) band (it corresponds to 3(T errors 
as taken from [10]) on the figure. However, as we can see on the left panel, for more of 
90% of the points the agreement is within 20%. This corresponds to the cases where the 
bottom-sbottom loop dominates as can be seen on the right panel, where we applied the 
cut sin(20^)Ai?Q^^^ > 0.02 to ensure that the bottom-sbottom loop is not negligible. 

So, in summary, we can say that the simplified approximate formulas of Eq. ( ^. and 
Eq. (^^) give a very good approximation of the full result for most of the cases. Only in 
the case when the bottom-sbottom loop is not the dominant diagram (small mixing) we 
get large deviations. 
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Figure 8: Comparison between the simplified analytical formulas and numerical results for the 
solar mass. In the left panel the LSP^^x" while in the right panel the LSP=f . The red (dark) band 
corresponds to the latest neutrino data [10]. 
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Figure 9: Comparison between the simplified analytical formulas and numerical results for the solar 
angle. In the left panel we present all points, and show in green the band that corresponds to more 
than 90% of the points. In the right panel we show all the points that have sin(26'jj)Ai?Q^^^ > 0.02. 
The red (dark) band corresponds to the latest neutrino data [10]. 



6. Rare radiative lepton decays 

As the parameters involved in the R-parity violating operator are constrained in order 
to predict neutrino masses in the sub-eV range, we have addressed [22] the question of 
whether this operator will induce rates for charged LFV processes of experimental interest. 
Some of them occur at tree-level such as double (3 decay [23, 24] and /j, — e conversion in 
nuclei [25]. One loop LFV decays as Ij li'j become interesting on this framework due to 
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the experimental interest in improving the current hmits [26]: 



BR{n 67) < 1.2 X 10 
BR{t < 1.1 X 10"^ 

BR{t 67) < 2.7 X 10"^ 



11 



(6.1) 



We have shown that the predictions for the last two processes are much lower than the 
above limits and will not constrain the BRpV model. For /i — > 67 the predictions are 
compatible with the current limit but could begin to constrain the model for the bounds 
that will be reached in current [27] or planned experiments [28], if only the atmospheric 
neutrino data were taken in account. However, the requirement that the one-loop induced 
AmgQL is in agreement with the solar neutrino data, implies that the predicted rates for 
/X — > 67 will not be visible [22], even in those new experiments. This can be seen in Fig. 10, 
where it is shown the contour plot for BR(/x — > 67) as well as the maximum values of ei, £2 
compatible with the solar neutrino mass. So, in conclusion, the experimental constraints 
on these rare decays do not constrain the BRpV model. 




0.001 



0.001 0.01 



0.1 1 
-ej(GeV) 



Figure 10: Contour plots for BR((U 67). The dashed line corresponds to m^^ = 0-01 eV. 



7. Probing Neutrino Mixing via SUSY Decays 

After having shown, in the previous sections, that the BRpV model produces an hierarchical 
mass spectrum for the neutrinos that can accommodate the present data for neutrino 
masses and mixings, we now turn to accelerator physics and will show how the neutrino 
properties can be probed by looking at the decays of supersymmetric particles. 

If R-parity is broken the LSP will decay. If the LSP decays then cosmological and 
astrophysical constraints on its nature no longer apply. Thus, within R-parity violating 
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SUSY, a priori any superparticle could be the LSP. In the constrained version of the MSSM 
(mSUGRA boundary conditions) one finds only two candidates for the LSP, namely the lightest 
neutralino or one of the right sleptons, in particular the right scalar tau. Therefore we will 
consider below, in detail, these two cases together with the possibility that a light scalar 
top can have sizeable R-parity violating decays. However, if we depart from the mSUGRA 
scenario, it has been shown recently [29] that the LSP can be of other type, like a squark, 
gluino, chargino or even a scalar neutrino. In section 7.4 we will briefly review their results. 



7.1 Probing Neutrino Mixing via Neutralino Decays 



If R-parity is broken, the neutralino is unstable and it will decay through the following 
channels: Xi ~^ T^ii^j^^k, i^iQQ, ^i^f^-k^ ^tlQ'^ ^il- It shown-*^ in Ref. [13], that 
the neutralino decays well inside the detectors and that the visible decay channels are 



quite large. This is shown in Fig. |11|, where in the left panel we have the neutralino 
decay length and on the right panel its invisible branching ratio. We see that the invisible 
branching ratio stays always below 10%, allowing most of the channels to be seen. We 
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Figure 11: In the left panel we show the neutralino decay length, hi the right panel it is shown 
the invisible neutralino branching ratio. ^ 

have shown in the previous sections that the ratios |Aj/Aj| and \ei/ej\ were very important 
in the choice of solutions for the neutrino mixing angles. What is exciting now, is that 
these ratios can be measured in accelerator experiments. In the left panel of Fig. |l^ we 
show the ratio of branching ratios for semileptonic neutralino decays into muons and taus: 
BR{x —>■ iiq'q)/BR{x — > Tq'q) as function of tan^^ATM- We can see that there is a strong 
correlation. The spread in this figure can in fact be explained by the fact that we do not 
know the SUSY parameters and are scanning over the allowed parameter space. This is 
illustrated in the right panel where we considered that SUSY was already discovered with 



^The relation of the neutrino parameters to the decays of the neutrahno has also been considered in 
Ref. [30]. 
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the following values for the parameters, 



M2 = 120 GeV, fi = 500 GeV, tan /? = 5, mo = 500 GeV, A = -500 GeV 



(7.1) 



We see that the correlatfon is now extremely good and a measure of those semileptonic 



10 



10 



10 



-2 




3, 



0.1 



510"^ 10"^ 



0.5 1 




tan^ e, 



aim 



Figure 12: Ratios of semileptonic branching ratios as a function of tan^ 6* a 
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On the left for 



random SUSY values and on the right for the SUSY point of Eq. (|7.l| ) 



branching ratios will be an important test for the model. By looking at other decays we 



can look at other neutrino parameters, as is indicated in Fig. 13. On the top row we have 
the correlation with U^^ and on the bottom row the correlation with the solar angle. As 
before, the left panels correspond to a random scan over the SUSY parameter space, while 
on the right panels we have assumed that SUSY was already discovered. As an example we 



took the SUSY parameters in Eq. (7.1) 



7.2 Probing Neutrino Mixing via Charged Lepton Decays 

After considering the case of the LSP being the neutralino, for completeness, we have also 
studied [15] the case where a charged scalar lepton, most probably the scalar tau, is the 
LSP. We have considered the production and decays of f, e and /i, and have shown that 
also for the case of charged sleptons as LSPs they will decay well inside the detector. This 
is shown in Fig. where on the left panel we show the production cross sections and on 
the right panel the decay lengths for all the sleptons. After showing that the sleptons will 
decay inside the detector we can correlate branching ratios with neutrino properties. This 
is shown in Fig. ^ where the branching ratios of the scalar taus are correlated with the 
solar angle, showing a strong correlation. 

7.3 Stop Decays and the Solar Angle 

It was shown [16] that the stop decays are complementary to decays as probes of 
the neutrino properties. The semileptonic decays of the stop can be important as it is 



shown in the left panel of Fig. 16, taken from Ref. [16], where it is plotted the tan/3 
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Figure 13: Correlation of neutralino decays with U^^ (top row) and with the solar angle (bottom 
row). On the left column we made a random scan over the SUSY parameter space, while on the 



right the parameters from Eq. (7.1) were assumed. 
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dependence of the branching ratio for the decay of ti into for several values of the 
neutrino mass. For m,^,^ = O.OGeV, the B(ti br) is still above 0.1% if tan/? is not 
too large. In the right panel of Fig. 16 we show the ratio of B(ti 6e^)/B(ti b i^i^) 
versus (ei/e2)^ for different values of cos^^. For definiteness we have fixed the heaviest 
neutrino mass at the best-fit value indicated by the atmospheric neutrino anomaly. One 
can see that the dependence is nearly linear even for rather small cos One sees from 
the figure that, as long as cos 9^ > 10~^ there is a good degree of correlation between the 
branching ratios into B{ti be~^) and B{ti — > b fi'^) and the ratio (ei/e2)^. Thus by 
measuring these branchings one will get information on the solar neutrino mixing, since 
tan^ ^soL is proportional to (£1/62)^ [11]- As a result in this model one can directly test 
the solution of the solar neutrino problem against the lighter stop decay pattern. This 
is also complementary to the case of neutralino decays considered in section 7.1. In that 
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tan^ 9q tan- 6q 

Figure 15: Ratios of branching ratios for scalar tau decays versus tan^ 9q. The left panel shown 
all data points, the right one refers only to data points with £2/63 restricted to the range [0.9,1.1]. 



case the sensitivity is mainly to atmospheric mixing, as opposed to solar mixing as can 



be seen from Figs. |12| and [IJ. Testing the solar mixing in neutralino decays at a collider 
experiment requires more detailed information on the complete spectrum to test the solar 
angle [13]. In contrast we have obtained here a rather neat connection of stop decays with 
the solar neutrino physics. 

7.4 Other LSP decays 

As we discussed in the introduction to this section, if we depart from the mSUGRA scenario, 
then the LSP can be of other type, like a squark, gluino, chargino or even a scalar neutrino, 
as it has been shown recently in Ref. [29]. As the decays of these LSP's will always depend 
on the parameters that violate R-parity and induce neutrino masses and mixings, it is 



possible to correlate branching ratios to the neutrino properties. This is shown in Fig. 17 



taken from Ref [29]. On the left panel we consider the case of chargino decays and show 
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Figure 16: On the left panel we plot the branching ratios for ii decays for mj:_^ = 220 GeV, 
/i = 500 GeV, M ^ 240 GeV, and = 100, 1 and 0.06 eV for cos 6*4 = -0.8. On the right panel 
we plot the ratios: B(fi fee+)/B(<i bfi^) as a function of (£1/62)^ for m^^ ~ 0.6 eV and 
I cos 6*4! > 0.1,0.01, 10-3. These plots were taken from Ref. [16] 



BR{x~^ — > ficc)/BR{x^ Tcc) as a function of the ratio (A2/A3)^. As we saw in section]^, 
this last ratio is correlated with the atmospheric angle. So, by looking at the chargino 
decays one can test the atmospheric mixing angle. In fact, as the authors of Ref [29] have 
shown, one can use the already very precise data on the atmospheric mixing angle to put 
bounds on the ratios of several branching ratios of chargino decays. On the right panel 
we show a very strong correlation for the solar mixing angle obtained [29] with the decays 
of squarks. Many other correlations can be obtained making the model over constrained. 
So, in summary, no matter what supersymmetric particle is the LSP, measurements of 
branching ratios at future accelerators will provide a definite test of the BRpV model as a 
viable model for explaining the neutrino properties. 



8. Conclusions 

The Bilinear R-Parity Violation Model is a simple extension of the MSSM that leads to a 
very rich phenomenology. Hopefully, it will be an effective model for the more theoretically 
attractive case where R-parity is spontaneously broken [31,32]. 

We have calculated the one-loop corrected masses and mixings for the neutrinos in a 
completely consistent way, including the RG equations and correctly minimizing the po- 
tential. We have shown that it is possible to get easily maximal mixing for the atmospheric 
neutrinos and large angle MSW, as it is preferred by the present neutrino data. We have 
also obtained approximate formulas for the solar mass and solar mixing angle, that we 
found to be very good, precisely in the region of parameters favored by this data. 

We emphasize that the LSP decays inside the detectors, thus leading to a very different 
phenomenology than the MSSM. In the mSUGRA scenarios, the LSP can be either the lightest 
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Figure 17: Left panel: ratios of branching ratios for chargino decays Dvdrsus (A2/A3)^; Right panel: 
ratios of branching ratios for squark decays versus tan^ ^SOL. These p^s were taken from Ref. [29]. 

neutralino, like in the MSSM, or a charged particle, must probably the lightest stau. In 
both cases we have shown that ratios of the branching ratios of the LSP can be correlated 
with the neutrino parameters. We have also shown, that the decays of the stop can be 
complementary to those of the LSP for the case of the solar mixing angle. 

In more general scenarios than the constrained mSUGRA, the LSP can essentially be any 
supcrsymmetric particle. However, also in these scenarios the neutrino properties can be 
correlated with ratios of branching ratios for these particles. 

So, in summary, no matter what supcrsymmetric particle is the LSP measurements of 
branching ratios at future accelerators will provide a definite test of the BRpV model as a 
viable model for explaining the neutrino properties. 
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